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(S Abstract 

In the framework of a specific scheme of the QCD sum rules for S- 
wave levels of the heavy quarkonium, one derives an expression, relat- 

[ — . ing the energetic density of quarkonium states and universal character- 

istics in the heavy quarkonium physics, such as the difference between 

f ^ ' the masses of a heavy quark Q and meson (Qq) and the number of 

heavy quarkonium levels below the threshold of (QQ) — ► (Qq) + (Qq) 
decay. 

Ph! Introduction 

43 

Powerful tools in studies of heavy quarkonia, the bound states of two heavy 

quarks, are phenomenological potential models [p], 0, §] and QCD sum rules 

[|] . An applicability of the approaches to the systems of two heavy quarks is 

caused by 

1) a low value of the ratio Aqcd/t^q ^ 1, where tuq is the heavy quark 
mass and A.qcd is the quark confinement scale, determining the inverse dis- 
tance between the quarks in the bound states, and 

2) a nonrelativistic motion of the heavy quarks inside the quarkonium, 

In the QCD sum rules, the low value of ratio Aqcd/wiq determines not 
large contribution of higher orders of the QCD perturbation theory over 
the quark-gluon coupling as ~ 1/ ^(tuq/Aqcd) <C 1 in the expansion of 
Wilson's coefficients, and it makes a suppression of nonperturbative quark- 
gluon condensate contribution, having the power form, as 0(Ag CD /mg) for 

1 



the gluonic condensate < as G > contribution into the sum rules for vector 



flU 

currents, for example. 

In the potential models, from data on the spectroscopy of the (cc) char- 
monium and the (bb) bottomonium one finds that the nonrelativistic quark 
motion allows one to get the phenomenological potential in the range of av- 
erage distances between the heavy quarks inside the quarkonia 

0.1 fm < r < 1 fm . (1) 

Being the potential of static sources for the gluon field, this potential 
must not depend on the flavours of sources. This flavour-independence is 
empirically confirmed for the QCD motivated potentials [|IJ. Such potentials, 
possessing different asymptotic properties in the regions of r — > and r — > oo, 
coincide with each other in the region ([I]), where they allow approximations, 
having a simple scaling behaviour. These approximations are the logarithmic 
[0] and power || laws 

(2) 
(3) 

By the virial theorem 

1 rlV 

(4) 
one finds 

< T L > = d L /2 = const. , (5) 

<T M > = j^i-CM + E), (6) 

where E is the binding energy of quarks in the quarkonium. Phenomenolog- 
ically, one has k <C 1, \E\ <C \cm\, so that in the region of average distances 
between the heavy quarks in the heavy quarkonium ([!]), the kinetic energy 
of quarks practically is a constant value, independent of the quark flavours, 

< T M >~ const. (7) 

Then from the Feynman-Hellmann theorem 

dE <T> 



V L (r) = 


c L + d L ln(A L r) 


i4/(r) = 


c M + d M (ALr) k 


<T >- 


1 dV 



du u 



where /i is the reduced mass of heavy quark system (QQ f ), one can get that 
the level difference in the system does not depend on the reduced mass of 
quarks, i.e. on the quark flavours, 

E(n, At) - E(n, //) = E(n, //) - E(n, //) . (9) 

Condition (^) means that the energetic density of heavy quarkonium levels 
does not depend on the quark flavours 

''" <P(n) , (10) 



dM n 



where <p(n) does not depend on /x. 

The described preprties of heavy quark potential are found phenomeno- 
logically. They cause the high accuracy of potential models for calculations of 
the heavy quarkonium masses with no account of spin-dependent splittings, 
5m{nL) ~ 30 MeV. 

The accuracy for the predictions of quarkonium wave functions in the 
framework of the potential models is low, for example, it is <5\I/(0)/\I/(0) ~ 
30 -7- 50%, since in this case the potential behaviour in the border points 
(r —>■ and r —>■ oo) becomes essential. 

In the QCD sum rules, the accuracy of predictions for the heavy quarko- 
nium masses is one order of magnitude lower than the accuracy of potential 
models, Stusr ~ 200 -r- 300 MeV. This fact is connected to that the consid- 
eration in the QCD sum rules takes a finite number of terms in the QCD 
perturbation theory for the Wilson's coefficients and a restricted set of the 
quark-gluon condensates, so that the results of such noncomplete considera- 
tion depend on an unphysical parameter, defining a scheme of the averaging 
in the QCD sum rules (the number of moment for a spectral density of current 
correlators or the Borel transform parameter). An additional uncertainty is 
related with a modelling of a nonresonant contribution into the current cor- 
relator, i.e. with the threshold of hadronic continuum. Such parametric 
dependences lead to the low accuracy of QCD sum rule predictions for the 
heavy quarkonium masses[|. Moreover, the use of weight functions, defining 
the averaging scheme and rapidly dropping with the energy rise, causes a 



1 The QCD sum rule accuracy in calculations of the leptonic constants (/,/,, /t) is 
higher (~ 20 4- 25%), since one uses the heavy quarkonium masses, known experimentally. 



suppression of the contribution of higher excitations in the quarkonium, so 
that, as a result, this contribution is neglected. 

Recently the QCD sum rule scheme has been offered in papers of refs. 
[|> S HI) so this scheme allows one to take into account the spectroscopic 
characteristics of higher S'-wave excitations. In these papers the following 
universal regularities have been derived. 

1) The scaling relation for the leptonic constants of S'-wave levels of the 
heavy quarkonium with the mass M and the reduced quark mass // is 

f - (^) 2 = const. , (11) 

that is in a good agreement (A/// ~ 5%) with the experimental data 
on the leptonic constants of T-, tp- and ^-particles, which are the quarkonia 
with the hidden flavours (QQ), so that Ajjl/M = 1, and one has || 

P 

— = const. , (12) 

independently of the heavy quark flavours in the (QQ) system. Relation (|P2]) 
essentially differs from the scaling law for the leptonic constants of heavy 
mesons (Qq), containing a single heavy quark, where in Heavy Quark Effec- 
tice Theory (HQET) [|l^] one has 

f ■ M = const. (13) 

Law (13) can be obtained from eq. ([TT|) in the limit jjl = m q m,Q / (m q + tuq) —>■ 
m q , niQ 3> m q , M — > tuq, so that \i does not depend on the heavy quark 
flavour. Eq.([U]) gives reasonable estimates for the leptonic constants of B- 
and D-mesons |Q, if one supposes \i ~ 330 MeV 0. 

2) The scaling relation for the leptonic constants of nS-levels in the 
quarkonium is 

f -k = -, (14) 

independently of the heavy quark flavours. Eq. fll4|) is in a good agreement 
with the experimental data on the leptonic constants in the families of i/j- 
and T-particles B (Af/f < 10%). 



3) The relation for the mass differences of nS-wave levels in the heavy 
quarkonium || is 

Mn-Mi Inn 

M^rwr^2' n - 2 > (15) 

independently of the flavours of heavy quarks in the quarkonium. Eq.flT5|) 
is in a good agreement with the experimental data on the masses of ip- and 
T-particles (5(AM)/AM < 10%), too. 

Regularities (Q), (P|), Ql5|) allow one to improve the accuracy of QCD 
sum rule results by one order of magnitude, approximately. This improve- 
ment takes place also for the heavy quarkonium masses, where the sum rule 
accuracy becomes comparable with the accuracy of potential models, so, in 
addition, the sum rules allow one to get the explicit relations. 

However, having derived relations (PD, (|14]), fll5|), one has used the 
phenomenological condition, stating the flavour-independence of the heavy 
quarkonium level density and coming from the analysis, made in the frame- 
work of the potential models. 

In the present paper, in the framework of the offered scheme of QCD 
sum rules, we derive the relation for the S-wave level density for the heavy 
quarkonium 

dn ' lnn^ 

where A = m^Qq) — tuq is the difference between the masses of heavy meson 
and heavy quark, nth is the number of S'-wave levels of the (QQ) heavy 
quarkonium below the threshold of quarkonium decay into the heavy meson 
pair (QQ) — > (Qq) + (Qq)- In the leading order, one has 

A = const. , (17) 

with the accuracy up to power corrections over the inverse mass of heavy 
quark |K|] (about the role of logarithmic and power corrections, see ref . ||1 1|| ) . 
In the leading approximation, stepping from the charmonium to the bot- 
tomonium, one can neglect a weak logarithmic variation of the number of 
levels below the thershold, 

In nth (66) - In n th (cc) . (18) 

From eq.(|16|)-(|18|) it follows that in the QCD sum rules one can show that 
the S'-wave level density of heavy quarkonium does not depend on the heavy 
quark flavours. 



Thus, the offered scheme of QCD sum rules allows one to do the quite 
complete consideration of heavy quarkonium and to use no external assump- 
tions, extracted from the phenomenological potential models, for example. 

In Section 1 we consider the scheme of QCD sum rules with the account 
of spectroscopic quantities for the heavy quarkonium and derive relation 
(16) . In Section 2 we make the phenomenological analysis of eq. (|T6|) . In the 



Conclusion the obtained results are summarized. 

1 Heavy Quarkonium Sum Rules 

Let us consider the two-point correlator functions of quark currents 

IV(g 2 ) = iJd 4 xe^ x <0\TJ^x)Jl(0)\0>, (19) 

U P {q 2 ) = i I ' d A xe iqx < 0|TJ 5 (x) 4(0)10 > , (20) 

where 

J M (z) = Qi(x) lfl Q 2 (x) , (21) 

Mx) = Qi(xhMx), (22) 

(23) 

Qi is the spinor field of the heavy quark with % = c, b. 
Further, write down 

IV = (-^ + ^)nv(? 2 ) + ^IW) , (24) 

where Hy and Us are the vector and scalar correlator functions, respectively. 
In what follows we will consider the vector and pseudoscalar correlators: 
U v (q 2 ) and U P (q 2 ). 

Define the leptonic constants fy and fp 

<0\J,(x)\V(\)> = ie^f v M v e lkx , (25) 

< 0\J 5ll (x)\P > = ik^fpe**, (26) 

where 

•V(aO = Qi(xh5ltiQ2(x) , (27) 



so that 

f A/f 2 
< 0\J 5 (x)\P >= % JP p e ikx , (28) 

mi + m2 

where |V > and |P > are the state vectors of l - and 0~ quarkonia, and 
A is the vector quarkonium polarization, k is 4-momentum of the meson, 

k 2 - M 2 

tvpy — IVlpy. 

Considering the charmonium (ip, ip' ...) and bottomonium (T, T', T" ...), 
one can easily show that the relation between the width of leptonic decay 
V — > e + e~ and fy has the form 

r<V - e+e~) = % e?aL &- , (29) 

where e^ is the electric charge of quark %. 

In the region of narrow nonoverlapping resonances, it follows from eqs. ( |T9| ) 
- (H) that 

-QmU^iq 2 ) = EfvnM 2 v J(q 2 -M Vn ), (30) 

71 

n 

i3mn£"V) = E/pnM^ n — ^-— ^% 2 -M^). (31) 
7r n (mi -t- m2j 

Thus, for the observed spectral function one has 

^m<f (g 2 ) = Ismn^tfj+p.ptf,^) , (32) 

TV 71 



where p(q 2 , /J. 2 ) is the continuum contribution, which is not equal to zero at 

> /x 2 - 

Moreover, the operator product expansion gives 



q 2 > /i 2 



n «^)(g 2 ) = n^ rt V) + C G {q 2 ) < ^G 2 > +C t (q 2 ) < mtQtQi >+..., 

71 

(33) 
where the perturbative contribution U^^' (q 2 ) is labeled, and the nonpertur- 
bative one is expressed in the form of sum of quark-gluon condensates with 
the Wilson's coefficients, which can be calculated in the QCD perturbative 
theory. 



In eq.(^) we have been restricted by the contribution of vacuum expec- 
tation values for the operators with dimension d = 4. For Cq (q 2 ) one has, 
for instance, |§] 

« P) = 1 f f3 (tf tl )(l-^ U.,%Htf + 3| 
G 192m 1 m 2 g 2 V 2w 5 1 - u w 4 / ' v ' 

where 

-2 2/ \2 2 4:171x1712 

q = q - (mi - m 2 ) , f = 1 ^ — . (35) 

The analogous formulae for other Wilson's coefficients can be found in Ref. [f|]. 
In what follows it will be clear that the explicit form of coefficients has no 
significant meaning for the present consideration. 

In the leading order of QCD perturbation theory it has been found for 
the imaginary part of correlator that || 

QmHy (q 2 ) = - — 2-(3ss — s 2 + 6mim 2 s — 2m\s)6(s — {m\ + rri2) 2 ($6) 

3mnf f) (g 2 ) = ^-(s-{m 1 -m 2 ) 2 )e{s-{m 1 +m 2 ) 2 ) 1 (37) 

where s = s — m\ + m|, s 2 = s 2 — Am^s. 

The one- loop contribution into Qmll(q 2 ) can be included into the consid- 
eration (see, for example, Ref.[U). However, we note that the more essential 
correction is that of summing a set over the powers of (a s /v), where v is 
defined in eq.(|35"D and is a relative quark velocity, and as is the QCD interac- 
tion constant. In Ref.j4| it has been shown that account of the Coulomb-like 
gluonic interaction between the quarks leads to the factor 

F <»> = T f 1-expMg) • (38) 

so that the expansion of the F(v) over as/v <C 1 restores, precisely, the 
one- loop 0(— ) correction 

„/ x 27r a* , 

F(v) « 1 - — -t . . . 39 

3 v 

In accordance with the dispersion relation one has the QCD sum rules, which 

state that, in average, it is true that, at least, at q 2 < 

i f^^y ds = 

n J s — q 2 



where the necessary subtractions are omitted. ^smli^ had ^ (q 2 ) and H^ CD \q 2 ) 
are defined by eqs.(^) - fl32|) and eqs.(^) - (|39"D , respectively, eg. (f40"D is the 



base to develop the sum rule approach in the forms of the correlator function 
moments and of the Borel transform analysis (see Ref.|4j]). The truncation 
of the set in the right hand side of eg. ([TOD leads to the mentioned unphysical 
dependence of the fpy values on the external parameter of the sum rule 
scheme. 

Further, let us use the conditions, simplifying the consideration due to 
the heavy quarkonium. 

1.1 Nonperturbative Contribution 

We assume that, in the limit of the very heavy quark mass, the power cor- 
rections of nonperturbative contribution are small. From eq. (|34]) one can see 
that, for example, 

4 P V) « 0(-}—) , A/m lj2 «l, (41) 

where v is fixed, q 2 ~ (m\ + rn.2) 2 , when 9mll' pert ' (q 2 ) ~ (mi + mq) 2 . It is 
evident that, due to the purely dimensional consideration, one can believe 
that the Wilson's coefficients tend to zero as 1/m 2 2 . 

Thus, the limit of very large heavy quark mass implies that one can 
neglect the quark-gluon condensate contribution. 

1.2 Nonrelativistic Quark Motion 

The nonrelativistic quark motion implies that, in the resonant region, one 
has, in accordance with eq.(|35"|). 



v -»• . (42) 

So, one can easily find that in the leading order 



3mnr tj (s) « %mll ( r rt \s) - ^- 2 s [^j , (43) 



so that with account of the Coulomb factor 



47r a s ,.., 

F(v) ~ , (44) 

O 1) 
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one obtaines 



^ rt) (^| S g) 2 . (45) 



1.3 "Smooth Average Value" Scheme of the Sum Rules 

As for the hadronic part of the correlator, one can write down for the narrow 
resonance contribution 



niTV) = /^E«*-M^), (46) 

J O U ,v> 



n «„, V) = I _±_ ? flnJ £k^ S(a - Mk) , (47) 



The integrals in eqs.(|46|)-(|47|) are simply calculated, and this procedure is 
generally used. 

In the presented scheme, let us introduce the function of state number 
n(s), so that 

n{m 2 k ) = k . (48) 

This definition seems to be reasonable in the resonant region. Then one has, 
for example, that 



-smU^ s \s) = sf Vn{s) Ts T,0(s- Ml k ) . (49) 

Further, it is evident that 

sp(-«J) = ^ !?»("-*)■ (50) 



and eq.(|4q) can be rewritten as 



The "smooth average value" scheme means that 

10 



It is evident that, in average, the first derivative of step-like function in the 
resonant region is equal to 

<^-j2 e ( n ~ k ) >- 1 • ( 53 ) 

k 

Thus, in the scheme one has 

<IirV)> - j^fvni/-^ 1 , (54) 

<IirV)> - / d \ s2f ^(s) dn^_ (55) 

J s — q 2 (mi + m 2 ) as 

Eqs.(^)-(^) give the average correlators for the vector and pseudoscalar 
mesons, therefore, due to eq.(|40"D we state that 



3m < rt hadr \q 2 ) >= QmU^ CD \q 2 ) , (56) 



that gives with account of eqs.(f45"D, fl54|) and (|55|) at the physical points 

Jl = as_ dMn f4/i\ 2 

where in the limit of heavy quarks we use, that for the resonances one has 

mi + m 2 w M , (58) 

so that 

fvn ^ fpn = fn ■ (59) 

Thus, one can conclude that for the heavy quarkonia the QCD sum rules 
give the identity of fp and fy values for the pseudoscalar and vector states. 
Eq.(^7]) differs from the ordinary sum rule scheme because it does not 
contain the parameters, which are external to QCD. The quantity dM n /dn 
is purely phenomenological. It defines the average mass difference between 
the nearest levels with the identical quantum numbers. 



Further, as it has been shown in ref.|T2[, in the region of average distances 



between the heavy quarks in the charmonium and the bottomonium, 

0.1 fm < r < 1 fm , (60) 
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the QCD-motivated potentials allow the approximation in the form of loga- 
rithmic law H with the simple scaling properties, so 

dn .„.. 

const. , (61) 



dM„ 



i.e. the density of heavy quarkonium states with the given quantum numbers 
do not depend on the heavy quark flavours. 

In ref . H it has been shown, that relation ( |B"TD is also practically valid for 
the heavy quark potential approximation by the power law (Martin potential) 
||], where, neglecting a low value of the binding energy for the quarks inside 
the quarkonium, one can again get eg. (|6"T|) . 

In ref.[g] it has been found, that relation (|6"TD is valid with the accuracy 
up to small logarithmic corrections over the reduced mass of quarkonium, if 
one makes the quantization of S"-wave states for the quarkonium with the 
Martin potential by the Bohr-Sommerfeld procedure. 

Moreover, with the accuracy up to the logarithmic corrections, as is the 
constant value. Thus, as it has been shown in refs.|5|, H], for the leptonic 
constants of S-w&ve quarkonia, the scaling relation takes place 

M U = m " St • < 62) 

independently of the heavy quark flavours. 

Taking into the account eqs.(|58D and fl59|) and integrating eqs. (|54D , (|55|) 
by parts, one can get with the accuracy up to border terms, that one has 

J dn dM n 7i \Mj V ; 



(64) 



Comparing eqs.(|57|) and (p3|), one finds 

dfn = 1 

f n dn 2n 

that gives, after the integration, 



&-%■ (65) 
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Table 1: The experimental values of leptonic constants (in MeV) for the 
quarkonia in comparison with the estimates of present model. 



quantity 


exp. 


present 


f<f> 

u 

/t 


232 ±5 
409 ± 13 
714 ± 14 


230 ± 25 
400 ± 40 
700 ± 70 



Table 2: The experimental values of leptonic constants (in MeV) for the 
nS'-bottomonia in comparison with the estimates of present model. 



quantity 


exp. 


present 


/i 


715 ±15 


input 


h 


487 ±16 


506 ± 10 


h 


429 ± 14 


412 ±8 


h 


320 ± 30 


358 ±7 


h 


369 ± 46 


320 ±7 


h 


240 ± 30 


292 ±6 



Relation (|65|) leads to that the border terms, which have been neglected in 
the writing of eq. (|63|) , are identically equal to zero. 

First, note that eq.(^), relating the leptonic constants of different quarko- 
nia, turns out to be certainly valid for the quarkonia with the hidden flavour 
(cc, 66), where 4/x/M = 1 @, § (see Table |XJ) . 

Second, eq. ( |62"D gives estimates of the leptonic constants for the heavy B 
and D mesons, so these estimates are in a good agreement with the values, 
obtained in the framework of other schemes of the QCD sum rules 0] . 

Third, taking a value of the lS'-level leptonic constant as the input one, 
we have calculated the leptonic constants of higher nS'-excitations in the 
charmonium and the bottomonium and found a good agreement with the 
experimental values pj (see Tables 0, |] and Figures [T], 0). 
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Further, from eqs. ([57]) and ( |65| ) it follows that 



f[_L = ^i 

n Mr, 7T 



4^y 

MJ 



dM n 
dn 



so that, neglecting the low value of quark binding energy (M n 
0(1/M))), one gets 



dM„ 



1 dMr, 



[n 



dn n dn 
Integrating eq.(^TJ), one partially finds eq.flTo]) 

M n - Mi In n 



1) 



M 2 -M x ' ln2 



and 



M 2 - M x 



dMn 
dn 



(n 



n> 2 



1) In 2 



(66) 
Mi(l + 

(67) 

(68) 
(69) 



Thus, in the offered scheme of QCD sum rules, one takes into account the 
Coulomb-like as/f-corrections and, neglecing the power corrections over the 
inverse heavy quark mass, one gets the universal relation for the differences 
of nS'-wave level masses of the heavy quarkonium. 

Eq. (^) for the differences of nS'-wave level masses of the heavy quarko- 
nium does not contain external parameters and it allows direct comparison 
with the experimental data on the masses of particles in the ip- and T-families 



Table 3: The experimental values of leptonic constants (in MeV) for the 
nS-charmonia in comparison with the estimates of present model. 



quantity 


exp. 


present 


h 


410 ± 14 


input 


h 


283 ± 14 


290 ± 10 


h 


205 ± 20 


237 ±8 


h 


180 ± 30 


205 ±7 


h 


145 ± 15 


183 ±6 
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Figure 1: The calculated dependence of nS - bottomonium leptonic constants 
and the experimental values of fr(nS)- 



Dependence (|68|) and the experimental values for the relations of heavy 
quarkonium masses are presented on Figure |3], where one neglects the spin- 
spin splittings. 

Note, the ^(3770) and t/>(4040) charmonium states suppose to be the re- 
sults of the 3D- and SS'-states mixing, so that the .D-wave dominates in the 
^(3770)-state, and the mixing of the 3D and 3S wave functions is accompa- 
nied by a small shifts of the masses, so that we have supposed M 3 = M^( 40 4o) . 

As one can see from the Figure ^|, relation (|68[) is in a good agreement 
with the experimental data. 

These facts show that the offered scheme can be quite reliably applied to 
the systems with the heavy quarks. 



Further, using eq.flBq), at q — one can write down 



s *h ds 



sfl 



dn 






! "i , d\nn 
ds 



ds 



fi ianth ■ 



(70) 



15 



fn, MeV 



400 








200 


1 


1 1 


• 
1 



1 2 3 4 5 n 



Figure 2: The calculated dependence of nS - charmonium leptonic constants 
and the experimental values of f^( n s)- 



From the other hand, in the leading approximation, one gets 



h og_ (M^\ 2 _ <^s_ (jp 



2tt KM) 2n \M 



[Sth ~ Si 



(71) 



and, further, 

sth ~Si~ 2MAE , (72) 

where AE = E t ^ — ttiq — ttiqi is the difference between the threshold energies 
for the decay (QQ 1 ) — > (Qq) + (Q'q) and the (QQ') pair production. 
In HQET p| one has 



AE = 2A + 0(1 /m Q ) 



(73) 



i.e. in the leading approximation one can take AE ~ 2A, being a constant 

value, independent of the heavy quark flavours. 

Then one finds 

/f = 05/4^ 2 JA__ (?4) 



M ix \M) In nth 



Comparing eq.(^) and eq.(|74]), one can easily find that in the leading ap- 
proximation 



dM n 
dn 



[n=l) 



2A 
In n t h 



(75) 
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a{n) 
3 




1 2 3 4 5 6 n 



Figure 3: The experimental values of nS - bottomonium (solid dots) and 
charmonium (empty boxes) mass differences a(n) = (M n — Mi)/(M 2 — Mi) 
and the dependence in the present model a(n) = Inn/ In 2. 



Having derived eq.(|75|), one has assumed, that 

1) the binding energy of quarks in the lS'-state is negligibly small, than 
the excitation energy of nS'-levels 



Ei < A ~ 



dM 



dn 



(76) 



2) the excitation energy of levels is small in comparison with the quark 
masses 



A 



dM 
dn 



< m Q 



(77) 



so that y^ ~ M, 

3) in the leading approximation the hadronic continuum threshold is de- 
termined by the messes of heavy mesons 

(78) 



s t h - M m) + M (Qlq) ~ m Q + m Q i + 2A , 
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4) the number of states below the threshold is finite and weakly depends 
on the heavy quark flavours. 

Thus, from eqs.(p^) and Q75D one can conclude that in the framework of 
the QCD sum rules, one gets the universal regularity for the density of S- 
wave quarkonium levels, so, in the leading approximation, the relation does 
not depend on the heavy quark flavours 

dM n , _. .„ . 

_ -JL( n = 1) = const. (79) 

an 

2 Analysis of Relation for Level Density 

Relation ( [7S|) is got in the leading approxomation over the inverse mass of 
heavy quarks, when one can neglect the spin-dependent splittings. Therefore, 
for the A estimate we will use the values of 5-level masses of the quarkonia 
(cc) &nd(bb). 

One can easily show 



m(n 3 Si) = m(nS) + ^Am{nS) , (80) 

3 

m(n 1 6*0) = m(nS) — -Am(nS) , 

where Am(nS) is proportional to the leptonic constant squared f% s |], ||, |[ 
i~2fl , so that from the previous Section it follows 



Am(nS) = ^1 . (81) 

From the experimental data one has Am^(lS) = 117 MeV, and, taking into 
account eqs. (|80| , |8l|) , one gets 

m f {lS) = 3.068 GeV , 

m^(2S) = 3.670 GeV , 

m r (lS) = 9.440 GeV , (82) 

m r (2S) = 10.012 GeV, 



where we have taken into account that 

A A as ^ 

asm 

and 



91n 2<T A f^ 

eff 



so that as(T)/as(ip) ^ 3/4 |13| . From eq.fl8~2|) one has 

(M 2 -M l % ~ 0.602 GW, (83) 

(M 2 -Mi)| T ~ 0.572 GeV , 

i.e. in average one has 

< M 2 - Mi >~ 0.587 ± 0.015 GeV . (84) 

In the T-family, where the leading approximation over the inverse heavy 
quark mass must be the most reliable, one has 

n th = 4 . (85) 

Then 

A =< M 2 - M x >~ 0.587 ± 0.015 GeV . (86) 

Estimate (^(J) of the important parameter in HQET is in a good agreement 
with the estimates, made in the QCD sum rules for the heavy mesons [TlJ 

A ~ 0.57 ±0.07 GeV . (87) 

However, one has to note that the offered estimate from the sum rules for the 
heavy quarkonium has the accuracy, that surpass the accuracy of estimate 
(|8~?D by about one order of magnitude and it is within the limits of the 
accuracy S A ~ 20 MeV, that can be achieved, because of the nonperturbative 



corrections in 


QCD [11|. 












From eq.(| 


(]) one finds the estimate 














(n = 1) ~ 0.85 ±0.02 
an 


GeV , 






(88) 


that is slight 1} 


' greater than the estimates, 


made 


: in papers 


of ref . 


M 


where 



dM/dn ~ 0.75 GeV was determined in the polinomial interpolation of heavy 
quarkonium masses. 
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Conclusion 

In the framework of the QCD sum rules, the expression for the density of 
S'-wave levels of the heavy quarkonium is derived 

dM n . 2A 

in = 1, 



dn ' In n t h 

that in the leading approximation does not depend on the heavy quark 
flavours. Here, A = ni(Qq\ —vtiq and n t h is the number of n^-levels below the 
threshold of (QQ') — ► (Qq) + (Q'q) decay. The analysis of the spectroscopic 
data on the charmonium and bottomonium allows one to do the estimate 

A ~ 0.587 ±0.015 GeV , 

that is in a good agreement with the recent estimates from the QCD sum 
rules for the heavy mesons, but the former has the better accuracy. 

The derived relation for the heavy quarkonium level density allows one to 
do a complete consideration of heavy quarkonium in the specific scheme of 
QCD sum rules with no use of the phenomenological data from the potential 
models. Such consideration essentially improves the accuracy of QCD sum 
rule predictions. 
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